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Q"^ ' We present results of the careful study of the Gribov copies problem in SU{2) 

^P . lattice gauge theory for the direct maximal center projection widely used in con- 

finement studies. Applying simulated annealing algorithm we demonstrate that this 
problem is more severe than it was thought before. The projected (gauge nonin- 
C^ ■ variant) string tension is not in the agreement with the physical string tension. We 

' , do not find any indications that P-vortices reproduce the full SU{2) string tension 

Q ' neither in the infinite volume limit nor in the continuum limit. 

> ■ 

: 1 Introduction 

H ; 

5^ , The idea that the center vortices are the objects responsible for confinement in the non- 

abelian gauge theories is rather old [|I|, |^. Recently it has been argued that the center 
projection might provide a powerful tool to investigate this idea |^. It is suggested that 
projection dependent P-vortices defined on the lattice plaquettes are able to locate thick 
gauge invariant center vortices and thus provide the essential evidence for the center vor- 
tex picture of confinement. So far 3 different center gauges have been used in practical 
computations: the indirect maximal center (IMC) gauge 0, the direct maximal center 
(DMC) gauge [^ and the Laplacian center gauge (see also for a new proposal). 

IMC and DMC gauges suffer from the Gribov copies problem |^ and this technical 
problem gave rise to several claims and counterclaims. At present we have a drama in 
five acts: 



i) Initially there was the claim [^, ^ that the projected string tension (the string tension 
which is due to P-vortices) reproduces the full SU{2) string tension. Thus P-vortices are 
responsible for the confinement of color in gluodynamics. 

ii) The problem of Gribov copies has been raised in Ref. ||^, where it has been 
demonstrated that there are gauge copies which produce P-vortices evidently with no 
center vortex finding ability since projected Wilson loops have no area law. At the same 
time these gauge copies correspond to higher maxima of the gauge fixing functional, than 
those used in Q]. 

Hi) In our previous publications P, [l^ we confirmed the existence of two classes of 
gauge copies: with nonzero and zero projected string tension az(2)- We resolved the 
problem raised in [^ since we found the copies with the highest maxima of the gauge 
fixing functional which correspond to nonzero crz(2) ■ At the same time the value of crz(2) 
was essentially lower than that obtained in Q] and in is disagreement with the physical 
string tension asu{2)- 



iv) In Refs. [|TT], Q it has been argued that the disagreement between az{2) and 
(^su{2) is due to strong finite volume effects and it disappears on large enough lattices. In 
Refs. [y, |1^, 0, the usual relaxation plus overrelaxation (RO) algorithm was used to 
fix the DMC gauge. 

v) Below we present results^, obtained with the help of more powerful gauge fixing 
algorithm, simulated annealing (SA), and show that the problem of low value of (Tz{2) 
persists even on large lattices with physical extension up to 3/?7i. We discuss here DMC 
gauge. Our results for IMC gauge will be presented elsewhere. 

2 Gauge fixing procedure 

2.1 Gribov copies 

Direct maximal center gauge ||^ in SU{2) lattice gauge theory is defined by the maxi- 
mization of the following functional: 



n,fj, ^ ' n,fj. 



with respect to local gauge transformations, Un,fi is the lattice gauge field, V is the lattice 
volume. Condition ([|) fixes the gauge up to Z{2) gauge transformation, and can be 
considered as the Landau gauge for adjoint representation. Any fixed configuration can 
be decomposed into Z(2) and coset parts: Un,^ = Zn^^Vn,^, where Z„^^ = signTiUn,^- The 

^The preliminary results are presented in Ref. pj] 



plaquettes Zn^^^ constructed from the links Z„^ have values ±1. The P-vortices (which 
form closed surfaces in 4D space) are made from the plaquettes, dual to plaquettes with 

Some evidence has been collected, that P-vortices in DMC gauge can serve to locate 
gauge invariant center vortices. It has been reported that the expectation value of 
the projected Wilson loops Wz{2), computed via linking number of the static quarks 
trajectories and P-vortices, have area law. The projected string tension az(2) is very close 
to the string tension of the nonabelian theory (Jsu{2)- This fact has been called center 
dominance. Another important observation was that the density of P-vortices scales as a 



physical quantity M, 13 



The main problem of DMC gauge fixing is that the functional F{U) (|l|) has many local 
maxima. This is the analogue of the Gribov problem in continuum gauge theories ^ and 
we call configurations corresponding to these local maxima Gribov copies. For gauge 
conditions with Gribov copies the usual expression for gauge dependent quantities WM 



does not provide unambiguous definition (see recent discussion of this subject in [|15|). It 
was shown that for some gauge conditions the gauge dependent quantities depend strongly 



on the local maxima picked up [0, [T^- Thus to remove this ambiguity it is necessary 
to find the global maximum or, if this is impossible, to approach the global maximum as 
close as possible]^. 

The usual algorithm for the gauge of the type (0) is the relaxation algorithm which 
performs maximization iteratively site by site. The relaxation becomes more effective if 
one uses the overrelaxation. It is also known that the SA algorithm is more effective than 
RO algorithm. It is very useful when ambiguities induced by Gribov copies become severe 



17| , |T9[. For this reason we use this algorithm for DMC gauge. 

We follow procedure proposed and checked in 0] : for given configuration we generate 
Ncop gauge equivalent copies applying random gauge transformations to the initial con- 
figuration, after that we fix the gauge for each gauge copy thus producing Ncop Gribov 
copies for each configuration. Then we compute a gauge dependent quantity X on the 
configuration corresponding to the Gribov copy with the highest local maximum of (|l]), 
Fmax{Ncop)- Averaging over statistically independent gauge field configurations and vary- 
ing Ncop we obtain the function < X{Ncop) > and extrapolate it to Ncop -^ oo limit. This 
should provide a good estimation for < X > computed on the global maximum unless 
the algorithm in use produces local maximum far from the global ones. 

Note that there exists another proposal [^ for the gauge fixing procedure which is 
free of gauge copies problem. In some particular limit this procedure corresponds to the 
search of the global maximum. There is also a class of gauge conditions ||2l|, |^, ^ which 
do not suffer from the gauge copies problem. 



^For another proposal to solve this problem see e.g. [n8| 



2.2 Simulated annealing algorithm 



The functional F(U) can be considered ||22| as the spin action: 



S{s) = F{W) = -^Y. '^'-^^ i^nUn,,S^a^(. + 1) (2) 

where the "spin variable" Sn, is the gauge transformation matrix in the adjoint represen- 
tation: Sn,ab = \tr{aagnCrbgi) ■ 

The lattice field Un,^i plays the role of (almost) random local couplings. Maximization 
of the functional F{U^) is equivalent to decreasing the auxiliary temperature T of the 
statistical system with the partition function, 

Z = J2^w(^S{s)) . (3) 

One starts with equilibrating this spin glass at high temperature. Subsequently the 
temperature, T, is decreased adiabatically. In the limit T — > the system approaches 
the ground state, i.e. the maximal value of S{s). Our SA implementation consists of the 
following three steps: 

1. Thermalization at T = 1.16. 

2. Gradual decreasing of T down to T = 0.02. 

3. Final maximization by means of the RO algorithm. 

In steps 1 and 2 an updating of the spin configuration was done with Metropo- 
lis algorithm. The temperature was lowered after every sweep by the quantity 6T{T). 
6T{T) has been tuned such that the spin action increased about linearly with the num- 
ber of iteration sweeps. This has been realized by subdividing full range [0.02,1.16] into 
38 intervals of width AT = 0.03. The corresponding differences of the average action 
AS(T) >= S(T) — S(T — AT) have been computed on equilibrated configurations and 
were found to be very stable against statistical fluctuations among different Monte Carlo 
configurations. We used the same 6T(T) for all volumes at given value of gauge coupling 
(5 and slightly modified it for different f3's. The number of sweeps to be performed within 
each interval (T — AT, T] was chosen to be proportional to AS(T) and, subsequently, 
the corresponding value of ST(T) has been determined. To approach close enough to the 
equilibrium we have to perform about 2000 sweeps at the step 2 (for step 1 only 20 sweeps 
was enough). Finally, at step 3 it is necessary to make O(IOO) sweeps of the RO algorithm 
to satisfy convergence criterion. In total this procedure was essentially more costly than 
the RO algorithm, also it gives the higher maximum of the functional (|I|) . 



3 Results 

3.1 Details of simulation 

Our computations have been performed on lattices 12^, IG'*, 20^ and 24^ for (3 = 2.4 and 
16^, 20^ for f3 = 2.5. The number of configurations was 150 for 12^, 50 for 16^, 30 for 20*^, 
20 for 24^ 

Using the described above gauge fixing procedure we calculate various observables as 
functions of the number of randomly generated gauge copies N^op- To make computations 
feasible we reduce the maximal value of N^op on large lattices, see Table 1. 

Table 1: The maximal values of Ncop for various lattices. 



/5 = 2.4 
/5 = 2.5 



L = 12 L = 16 L = 20 L = 24 



10 7 7 4 

- 7 7 7 



The projected Wilson loops Wz{2)(yC) are defined as 

Wz(2){C) =< exp{27r£(Ep,C)} > 

where £(Sp, C) is the 4D linking number of the closed surface, Sp, formed by P-vortex 
and closed loop C. We used two ways to estimate projected string tension (7z(2)'- from 
the Creutz ratios 

,, n 1 M^z(2)(/.-f)W-z(2)(/+l,/+l) ,,, 



and fitting static potential 



Vz{2){R) = - lim log 



Wzi2)iR,T + l] 



T^^ " Wzi2){R,T) 

to the usual form 

Vo + c/R + azi2)a^R. (5) 

The P-vortex density p is another quantity of interest. It is defined as follows 

1 



3.2 Maximized functional and extrapolation to N, 



cop 



OO 



In DMC gauge the average value < Fmax > is an important indicator of the quahty of the 
gauge fixing since gauge noninvariant observables we are measuring are strongly correlated 
with < Fmax >■ In Fig.l we plot < Fmax > as a function of Ncop- For comparison we 
show the results for RO algorithm and the algorithm suggested in ^j : before applying RO 
algorithm the random gauge copy was fixed to Landau gauge. We call it LRO algorithm. 
One can see that SA algorithm indeed produces highest maxima. We checked that this is 
true for other lattices as well. 
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Figure 1: The average value of the highest local maxima < Fmax > for 3 algorithms as 



a function of A^, 



copi 



16, /5 = 2.5. The solid lines are the fits using the asymptotic 



expression ( [A. 1| ) 



One can draw two conclusions from Fig.l. First, it is clear that the local maxima 
generated with LRO procedure of |§] can be safely ignored. Second, RO algorithm does 
not permit to reach the global maxima and results obtained with this algorithm must be 
revised with the help of SA algorithm. As we mentioned the gauge dependent quantities 
are strongly correlated with the value of Fmax- 

The important question is how close we approach the global maximum. Theoretically 
SA algorithm permits to find the global maximum. In practice this is complicated task 
because of large autocorrelation length at small temperature T. We checked for a few 
configurations that increasing the number of sweeps in our SA algorithm by factor 5 



does not lead to any essential change in F^ax- This indicates that we are close to the 
equilibrium during SA cooling and thus our results are close to the global maximum when 



we make the extrapolation to N^ 



cop 



cxo limit. 



As we show in the Appendix there are 2 possibilities for dependence of < Fmax > on 
Ncop- Our data show that at small Ncop eq. (|A.2|) provides better fit while at large Ncop 



eq. (|A.l|) is better. For this reason we use eq.( |A.l| ) to extrapolate the data to N^op -^ oo 



limit. In Fig.2 we depict Xz(2)(5,5) and < p > as functions of < F^ax >■ The linear 
dependence seen in this figure proves that the use of the same fitting function (eq. ( |A.1| )) 
for < Fmax > and for Xz{2) {I, I) and < p > is indeed well grounded. 
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Figure 2: Xz(2)(5, 5) and < p > vs. < F^ax > for L = 20, /? = 2.4. 



3.3 P-vortex density 

The physical density of P-vortices < p/a^ > is the important indicator of their relevance 
for the infrared physics. If, e.g., it goes to zero in the continuum limit, then P-vortices 
cannot play any important role. It is known that for the completely uncorrelated vortices 
< p/a^ >= i(T^(2) PI- 

To compare the scaling of density with that of the unprojected string tension we show 



in Fig. 3 the ratioQ < p > /o-su{2)<^^- One can see that for two [3 values this ratio is 
almost the same, i.e. scaling properties of < p/a? > and asu{2) are similar. The volume 
dependence is rather weak. Our value for < p > is essentially lower than that obtained 
in |jl2| and this difference increases with f3. It will be clear from our results in the next 
section that < p/a"^ >~ crz{2), i-e. P- vortices are far from being uncorrelated. On the 
other hand (at least on the lattices up to 24^) the density of P- vortices does not reproduce 
the full string tension, (Tsu(2), as it was suggested in Ref. [^. 
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Figure 3: The ratio < p > /crsu{2)Ci^ vs. Ncop- The solid lines are the fits obtained using 



asymptotic expression (|A.l 



3.4 Projected string tension az{2) 

3.4.1 Creutz ratio Xz{2) 

At first we discuss the final volume effects. In Fig. 4 we show the dependence of Xz{2) (4, 4) 
on Ncop for jS = 2.4. From the figure one can see that finite volume effects are essential 
for L = 12 only. For the other values of / we obtained the similar behavior. At /? = 2.5 
we found a sizable finite volume effects only on L = 16 lattice. 



^CTs(7(2)a^ 



0.0728(6) for /3 = 2.4|3|, 0.0350(4) for /3 = 2.5 |24| 
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Figure 4: The dependence of the Creutz ratios x^(2)(4, 4) on the number of gauge copies 
Ncop for P = 2.4. The sohd hues are fits using eq. ( |A.1| ). 



Usually one uses the Creutz ratio Xz{2){I, I) (i|) to estimate the projected string tension 
because the dependence on I was found to be rather weak. In our study we found that 
it was not quite true. In Fig. 5 we plot the ratio xz{2){I,I)/o'su{2)a'^ as a function of the 
dimensional distance ^r = I ■ a. The data shown were obtained for maximal values of 
Ncop given in Table 1. The values extrapolated to Ncop — > C)0 limit are only slightly lower. 
One can see that Xz{2) {I, I) is not constant but increases with I at least in the range 
0.3/m ^ I ■ a ^ 0.7fm. The similar behavior can be seen in Fig.3 of |T^. We found 



that this effect weakens, when Ncop increases, but does not disappear completely. It is 
evident from this figure that even on our largest lattice (L = 24 at /5 = 2.4), which has 
extension in physical units of about 3/m, Xz(2){I, I) is much lower than the unprojected 
string tension. This invalidates the suggestion of |jl2| that Xz{2){I, I) approaches <ysu{2)0''^ 
on lattices larger than the size of Z(2)-vortex (~ Ifm). Let us remind that the authors 
of Ref. |]12[ have used the RO algorithm to fix the gauge, while we are using the SA 



algorithm. 



*The physical distance scale is set by the value (Jsu{2) — (440 Mev^ 
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Figure 5: The dependence of Xz(2){I,I)/o'su{2)(i'^ on the distance in physical units. N, 
are the largest values used for the given lattice (see Table 1). 
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Our estimation for o-z{2)/o'sui2) made from projected Creutz ratio measured at dis- 
tances 0.7 — 0.9fm is 0.66(2). Comparing our data obtained at given f3 on lattices of 
various size (see Fig. 5), we can conclude that for lattices of physical extension larger than 
1.7 fm there are no sizable finite volume effects for the string tension extracted from the 
Wilson loops of the size < Ifm. Still we cannot exclude that the projected string tension, 
measured from projected Wilson loops with a very large extension, approaches the full 
string tension. 



3.4.2 Static potential \/z(2)(i?,T) 

In view of the dependence of Xz{2){I, I) on I it is interesting to study the projected static 
potential Vz(2){R) = ^i'^T^ooVz(2){R,T), Vz(2){R,T) = -log w^I^(r.t) ■ ^^ i^ i*^ ^^^"^ 



Z(2) 



(R,T) 



from Fig. 6 Vz{2){R,T) rises as a function of T. If we assume that for Wz{2){R,T) the 
usual expansion 

Wz(2){R,T) = Co{R)exp{-Vzi2){R)T) + C,{R)exp{-Vi{R)T) + ... (6) 

is valid, then such behavior implies that Ci{R) < 0. 

To compute Vz(2) (R) we first extrapolated Vz(2) {R, T) as a function of N^op to Ncop 
limit and then evaluated Vz{2){R) using large T values. The results are depicted in Fig. 7. 
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Figure 6: Vz(2) (-R, 7") as a function of T for several values of i?. L = 24, /3 = 2.4, A^^ 
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Figure 7: The static potential Vz(2){R)- L = 24, /5 = 2.4. The solid line is fit using 
eq.(l). 
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The value of (Jz{2) obtained from the projected static potential is in agreement with 
that obtained from the projected Creutz ratio (Section 3.4.1). 

4 Conclusions 

We made the most careful up to date gauge fixing of DMC gauge employing SA algorithm. 
We have shown that this algorithm permits to get the higher local maximum of the 
functional F{U) than that obtained before with other algorithms. Though increment in 
< Fmax > looks uot impressivc it gives rise to substantial change in the values of gauge 
noninvariant observables. This, in turn, brought us to conclusions different from those 
made in [|T2|. We computed the projected string tension az(2) and P- vortex density < p > 



on volumes with physical extension from 1.4/m up to 3/m. Our results imply that only 
on our smallest lattices {L = 12 for (3 = 2.4 and L = 16 for f3 = 2.5) there are moderate 
finite volume effects. Our data for lattice size larger than 1.7 fm show the absence of the 
finite volume effects. 

The ratio of P- vortex density to <Jsui2) weakly depends on j3 (see Fig. 3) and the early 
asymptotic scaling of this density (observed in [§|, |12|) does not exist. 

The projected string tension, az{2), does not reproduce the full string tension, asu{2)- 
The obtained value of the ratio crz(2)/o'sui2) ~ 0.66(2) is rather low. Comparing results 
for two lattice spacings we do not see any convergence of this ratio to 1 in the continuum 
limit (further studies at larger values oi jS on larger lattices are important). Let us remind 
that abelian monopoles are responsible for 92% — 94% of the full SU{2) string tension 



17, 25 



Since the P-vortices in DMC gauge do not reproduce the full SU{2) string tension 
they cannot be responsible for confinement. On the other hand the fact that az{2) is of 
the same order as <7su{2) probably means that the DMC gauge is close to some "ideal 
gauge" which detect perfectly B the thick gauge independent center vortices EH]. Recent 



results [^ show that an example of such gauge is the Laplacian Center Gauge which is 



free of Gribov copies problem and az{2) ~ o'su{2) in this gauge. 
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Appendix A 

Let the probability distribution of the random variable x is p{x) and let p{x) is nonzero 
on the interval bo — 6i,xo + 62]- Let us introduce 



P{x) = dy p{y) 



'xo~5i 

with properties 

P'(x)=p(x), Pixo-5,) = 0, Pixo + 62) = l. 

Then probability distribution P^'^\x) to get value x as maximal out of N attempts is 
equal to 

p(i)(a;) =p{x) 

P^^\x) = CNPixf-^x) 

where normalization constant 

C],' = / dx P{xf~'p{x) = / dP P^-' = j-PX = ^ 

We have to evaluate 

rxo+52 rxo+52 

x(^) = / dx xP^^\x) =N dx xp{x)P^-\x) 

J xa^&i J X{)—5i 

/&2 
dy yp{y)P^'^{y) {x = y + xq) 
-<5i 



'' ...JP^'iy) _^ , ..o^,„.„fe ['' ,..T.N, 



= Xo+ dy y—r^ = Xo + yP^ (y) |\ - / dy P'' (y) 

= Xo + 62- r dyP''{y) = x- r dy P^ {y) 

J-5i J-5i 

where x = xq + 62 is evidently the value of x^^' in the limit A^ ^ 00. The integral in the 
last equation can be evaluated for large N in the following way: 

13 



/ = -/'' dy P'^iy) = - f' dy e^'°3(^(2')) 
We expand the function f{y) = log{P{y)) at its maximum y = 62'- 

Then 

I = - [ de e^(P(^2)e+ib'(52)~p2(52)]e2 + ...) 

Let us consider 2 cases: 

1. P(<52) 7^ 0. 



1 / . «^A_^,_L J_^„','A„^_^2/r„^u2_| 1 /, „ / 1 



Thus in this case the fitting function for numerical data should be: 



/i(A^) = ^ + § + ^ + - (A.l) 



2. p{62) = 0. 



de e3[p'('52)-p'('52)]^ 



1 / . i^„'(-A„^_„2/r„^u2_ 



Then fitting function is 



N J-^(5i+52) 



MN)=x + ^ + ^ + ... (A.2) 



If ^((^2) is a very small number then there might exist interval of A^ where eq. ([A.2|) 
provides better fit, while in the limit of large N eq.( [A.l| ) is valid. The data for our 



observables (< F^ax >, < P >, Xz(2)iI,I)) indicates just this situation. For this reason 
we made fits with eq.( |A.lD discarding small values of A^^ 



'cop- 
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